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Abstract 

We investigated the transverse beam dynamics in a thermal wave model by using 
a functional method. It can describe the beam optical elements separately with a 
kernel for a component. The method can be applied to general quadrupole magnets 
beyond a thin lens approximation as well as drift spaces. We found that the model 
can successfully describe the PARMILA simulation result through an FODO lattice 
structure for the Gaussian input beam without space charge effects. 
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The thermal wave model is an efficient way to study the beam dynamics of relativistic 
charged particles. The Schrodinger-type equation in the model governs the time evolution 
of the beam wave function whose squared magnitude is proportional to the particle number 
densitiespTj. The model has successfully explained the filamentation of a particle beam and 
the self-pinching equilibrium in collisionless plasma [2 J. It was also used to estimate the 
luminosity in a linear collider where a spherical aberration was present [3]. The model can 
also provide some insight into a halo formation by introducing a Gaussian slit [4]. 

Transverse beam dynamics in a one spatial dimension is another application area of the 
thermal wave model. In Ref . [5] , the authors investigated the beam wave function through 
a quadrupole magnet with sextupole and octupole perturbations followed by a long drift 
space under a thin lens approximation. There is also a paper on the phase space behavior 
of particle beams in the transverse directions where the Wigner and Husimi functions are 
used as the phase space distribution functions^. 

In this work, we investigate the transverse beam dynamics in a two dimensional trace 
(x — x' or y — y') space in the thermal wave model by using the functional integral method 
[7] . Because the method can be extended to general lattice structures including quadrupole 
magnets beyond a thin lens approximation limit and it can treat the beam optical ele- 
ments individually, it is possible to systematically analyze a beam motion in a realistic 
environment such as an FODO lattice. We found that the model can successfully explain 
the PARMILA[8] simulation results with Gaussian input beams in a two dimensional trace 
space under the condition that the space charge effects are negligible. We note that this 
method can explain a low energy particle behavior as well as the relativistic motion of 
the charged particles if the important interactions are related to the external linear optical 
elements such as the quadrupole magnets and the random motion described by a beam 
emittance. 

In the thermal wave model, the time evolution of the beam wave function for the rela- 
tivistic charged particles can be described by the Schrodinger-type equation in the trans- 
verse directions. Because the beam dynamics are usually described in the two dimensional 
trace (x — x' or y — y') space, it is important to see whether the one dimensional version 
of the equation can explain the beam dynamics in the projected space or not. The one 
dimensional Schrodinger-type equation in x direction is given by 

where z = ct is the longitudinal distance of the beam movement and U x (x, z) = u x (x, z)/mo r ) r 0^c 2 
is the dimensionless potential with the relativistic parameters, (3 T = v/c and 7 r = (1 — 
Pr) -1 ^ 2 - The parameter e x is related to the emittance of the particle distribution in the 
x x space, which is explained later. The transverse particle distribution can be obtained 
by the squared magnitude of the beam wave function, N \ip x {x, z)\ 2 with the particle num- 
ber of N. In this convention, the beam wave function satisfies the normalization condition 
as follows, /f^ \ip x {x, z)\ 2 dx = 1. The corresponding equation for the time evolution of the 
wave function in the y direction can be obtained by replacing x with y in Eq. (pQ). We note 
that the time evolution of the beam wave function in y direction is independent of that in 
the x direction because the one dimensional Schrodinger-type equation in the y direction 
includes the different parameter of e y and we considered the linear external forces only. In 
the following analysis, the parameters and the functions in the x direction are used without 
subscript if there is no confusion. 

We can solve the differential equation by imposing the following two boundary condi- 
tions, a 2 (z = 0) = dp and = [5]. The a denotes the root mean square (rms) size 
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of the beam distribution and p is the curvature radius of the beam wave function along the 
beam direction. 

Another efficient way to solve the differential equation is known as the functional inte- 
gral method [7J where the resulting wave function is given by the product of a kernel (or 
propagator) and the initial beam wave function, 

/oo 
dxiK(x f , z f ; x h Zi) ^{x h Zi). (2) 
-oo 

Since a kernel represents an optical element like a quadrupole magnet or a drift space, the 
functional method can separate a mult i- components problem into several single-component 
problems. This property is the main advantage of this functional method in the thermal 
wave model. 

We can obtain the kernels from the path integral method [7] directly as follows, 

K(x f , z f ; Xi, Zi ) = J V[x(z)\ e iS ^, (3) 

where S(z) = f* f dzC(x(z),x'(z)) is called the action. The Lagrangian, C, of a system is 
the difference of the kinetic and potential energy terms. 

In this work, we will restrict our attention to a system consisting of quadrupole magnets 
and drift spaces. The potential energy terms of the beam optical elements are given by 

jj, x _ J for a drift space, 

1 \k\x 2 for a focusing quadrupole magnet, 

where k\ is positive in the focusing case. The potential term for the defocusing magnet is 
-(A:i/2)x 2 . 

The kernel, K , for a drift space which has no potential term is given by 

«**„v,«,«) = (^^)) 1/2 ^ ( "" I,)2 - W 

The kernel, Kf, for the focusing quadrupole magnet is given by 

/ n— \ 1/2 _ 

i / \ I V "'l \ i \ 1 \ (x 2 t +x 2 ) cot \Zk\z— 2x tX4 esc \/kiz\ /p \ 

K f(xf,Zf',Xi,Zi) = ■(/!-( e 1 f J ' ( 6 ) 

\27ue sm[\fki[Zf — Zi)) J 

For the defocusing case, the kernel is obtained easily by replacing the cot and esc functions 
in Eq.([6]) with coth and csch functions, respectively. 

Since the potential energy terms are related to linear forces only, the integration in Eq. 
(j2J) becomes very simple if the initial beam wave function is a Gaussian-type such as 
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where o\, p\, 9\ are the initial values of the rms beam size, the curvature radius, and the 
input phase, respectively. 

After the input beam passes through a linear optical element, the beam wave function 
remains the Gaussian-type such as 



ifa(x,z) = ( - \. ) exp 



x 2 .fx 2 



Aa 2 {z) \2ep 2 (z) 



6o(z) 
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The different forms of the parameter functions, $2(2), 02(2), and p2{z), characterize the 
properties of each optical element. 

In a drift space, the functions are given by 
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In a focusing quadrupole magnet, they are given by 
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with o"o = e/ (2y^i). For a defocusing lens, the functions can be obtained by replacing \fk 



for the focusing case with iy/ki. We can easily check to see if Eq. (jHJ) is the solution of 
Eq. ([I]) by inserting the obtained beam wave function into the differential equation. 

First of all, we studied how to relate the model parameters, cri, pi, e, of the input Gaus- 
sian wave function in Eq. ([7]) to the twiss parameters and the unnormalized rms emittance, 
«i, Pi, e rms . Because P is defined as <r 2 /e rms for the Gaussian distribution, we obtained 
o"i = V e rmsPi- From the definition of 1/p = (l/a)(da/dz) [5], we can easily obtain 
1/pi = — aij Pi where we used dp/dz = —2a. Motivated by the quantum mechanical 
relation between the wave functions in the configuration and momentum spaces, we defined 
the wave function in the x' space as the Fourier transformation of the Gaussian beam wave 
function as follows, 
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with tan(20id) = 2crf/(epi). The initial particle distribution in the x' space is proportional 
to |0i(x')| 2 . Because we can define 7 as crj/e rms for a Gaussian distribution in x' space, 
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we obtain a\a\ d = e 2 rms (3i^i. Comparing it with a\o\ d = e 2 /4 + 

trms&i which can be 

obtained from Eq. ( fTBT) . we can obtain e = 2e rms where we used the relation between twiss 
parameters, fliji — a\ = 1. We also obtained 1/pid = — ai/71 from Eq. (JTTj). We note that 
the relations between the model and physical parameters are valid for the wave functions 
at each of the beam optical elements. 

We note that the above analysis for the time evolution of the beam in x — x' space is 
also valid in the y — y' space if we use the one-dimensional Schrodinger-type equation in y 
direction with the emittance parameter of e y . In the following analysis, we studied the time 
evolution of the beam wave functions in both the horizontal (x — x') and vertical (y — y') 
spaces by using the equations with different emittance parameters, e x and e y . 

In order to check on the validity of the solutions, we compared them with the PARMILA 
simulation results with 50,000 macro particles through the FODO lattice in the horizontal 
direction. It corresponds to the DOFO lattice in the vertical direction. The field gradient 
and effective length of the quadrupole magnets in the lattice are 10.0 T/m and 0.2 m, 
respectively. The length of the drift spaces is 0.5 m. The particle type is proton with a 
kinetic energy of 100 MeV. We selected a random distribution of the particles in the trace 
spaces and neglected the space charge effects. The weighting function of the distribution 



is a Gaussian-type truncated at four times the standard deviation. Figure 1(a) (2(a)) 



and Figure 1(b) (2(b)) show the particle distributions of the input beam in the x(y) and 
x'(y') directions of the horizontal (vertical) trace space, respectively. The histograms are 
the PARMILA result with 50,000 macro particles. The real lines represent the Gaussian 
input beam for the model calculation. They are obtained by fitting the histograms of the 
PARMILA results. In the all figures of this work, we used the same normalization factors 
of the distribution functions as those of the input functions. We found that the beam wave 
function of Eq. (|T5l) describes the initial particle distribution very well in both x' and y' 
directions. 

The properties of the input beam are summarized in Table [TJ From the relations between 
the model and physical parameters, we can obtain the input values of the model parameters 
as follows, 



o\ x = 0.76 mm p\ x = —1.48 m for the horizontal direction 
o\ y = 0.55 mm p\ y = —0.40 m for the vertical direction 
ki = 6.74 m -2 , 



where ki = qG / {m^ r [3 r c) with the quadrupole field gradient, G. 

Figure [3] and Figure H] show the particle distributions at the end of each optical element. 
The histograms and real lines represent the PARMILA simulation results and the model 
calculations in the x and x' directions, respectively. Since the beam wave functions at 
each stage are Gaussian-type in the x direction, the wave functions in the x' direction can 
be obtained by applying Eq. (|T5l) . Corresponding figures in the y — y' space are given in 
Figure [5] and Figure [6J We note that the distribution functions are proportional to the 
\^ix(x)\ 2 (\i>i y (y) | 2 ) and \<p lx (x')\ 2 (\4>i y {y')\ 2 ) in x (y) and x' (y') directions, respectively. 
From Figure 0] and Figure El we can conclude that the Fourier transformation is a valid 
method to define the wave functions in the divergence directions. The figures show that the 
one dimensional Schrodinger-type equation of thermal wave model describes the PARMILA 
simulation result through the FODO (or DOFO) lattice successfully in the two dimensional 
trace (x — x' or y — y') space. In order to check on the result quantitatively, we compared 
the rms beam sizes obtained by the model with the values obtained by the best-fit of the 
PARMILA result. It is summarized in Table [2j It shows that the model results are the 
same as the simulation to within 0.8 %. 
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In conclusion, we studied the one-dimensional Schrodinger-type equation in the thermal 
wave model which describes the beam behavior in the x — x' or y — y' spaces. Some 
relations were obtained between the model parameters and physical parameters such as the 
twiss parameters and unnormalized rms emittance. We used a functional method in order 
to solve the differential equation with the Gaussian input distribution under the condition 
of the negligible space charge effects. The main advantage of this functional method is 
that we can calculate the effects of each beam optical element separately. The information 
of each element is summarized in a kernel. The final beam wave function of one optical 
element is obtained easily by the Gaussian integration of the product between the kernel 
and the initial Gaussian wave function. We found that there is a good agreement between 
the PARMILA simulation and the model calculation if we neglect the space charge effects. 
Even though there are some limits to the application of this method, this functional method 
is a very efficient tool to study the transverse beam dynamics in the thermal wave model. 
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Table 1: The twiss parameters and unnormalized rms emittances of the input beam in the 
horizontal and vertical directions. 





a 


j3 (m/rad) 


e rms (10~ 7 m-rad) 


horizontal (x) axis 


1.62 


2.41 


2.46 


vertical (y) axis 


2.87 


1.16 


2.57 



Table 2: The rms beam sizes obtained by the model and the best fit of the PARMILA 
result in the horizontal and vertical (numbers in parentheses) directions. 





model (mm) 


PARMILA (mm) 


deviation (%) 


After a F(D) lattice 
After a drift space 
After a D(F) lattice 
After a drift space 


0.572(0.353) 
0.251(0.369) 
0.554(0.460) 
1.521(0.644) 


0.575(0.355) 
0.253(0.370) 
0.558(0.461) 
1.530(0.647) 


-0.52(-0.56) 
-0.79(-0.27) 
-0.72(-0.22) 
-0.59(-0.46) 
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(a) 



(b) 



Figure 1: Particle distributions of the input beam in (a) x-direction and (b) x'-direction of 
the horizontal trace space. The histograms and real lines represent the PARMILA results 
and the model predictions, respectively. 
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divergence angle in y-direction (mrad) 



(a) 



(b) 



Figure 2: Particle distributions of the input beam in (a) y-direction and (b) y'-direction of 
the vertical trace space. The histograms and real lines represent the PARMILA results and 
the model predictions, respectively. 




Length in x-direction (cm) Length in x-direction (cm) 

(c) (d) 

Figure 3: Particle distributions on the horizontal axis with histograms for the PARMILA 
results and real lines for the model calculations: (a) after a focusing quadrupole (b) af- 
ter a drift space (c) after a defocusing quadrupole (d) after a drift space using the same 
normalization as the input beam distribution function. 
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divergence angle in x-direction (mrad) divergence angle in x-direction (mrad) 

(c) (d) 

Figure 4: Particle distributions on the horizontal divergence axis with histograms for the 
PARMILA results and real lines for the model calculations using Fourier transformation of 
the beam wave function : (a) after a focusing quadrupole (b) after a drift space (c) after 
a defocusing quadrupole (d) after a drift space using the same normalization as the input 
beam distribution function. 
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Length in y-direction (cm) Length in y-direction (cm) 

(c) (d) 

Figure 5: Particle distributions on the vertical axis with histograms for the PARMILA 
results and real lines for the model calculations: (a) after a defocusing quadrupole (b) 
after a drift space (c) after a focusing quadrupole (d) after a drift space using the same 
normalization as the input beam distribution function. 
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divergence angle in y-direction (mrad) divergence angle in y-direction (mrad) 

(c) (d) 

Figure 6: Particle distributions on the vertical divergence axis with histograms for the 
PARMILA results and real lines for the model calculations using Fourier transformation of 
the beam wave function : (a) after a defocusing quadrupole (b) after a drift space (c) after 
a focusing quadrupole (d) after a drift space using the same normalization as the input 
beam distribution function. 
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